Wigner function of pulsed fields by direct detection 
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q-i We present the reconstruction of the Wigner function of some classical 
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pulsed optical states obtained by direct measurement of the detected-photon 
probability distributions of the state displaced by a coherent field. We use 



a photodetector endowed with internal gain, which is operated in the non- 
■ photon-resolving regime. The measurements are performed up to mesoscopic 

intensities (up to more than 30 photons per pulse). The method can be 
O applied to characterize nonclassical continuous-variable states. © 2009 
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The effective implementation of quantum communication protocols requires the generation 
and measurement of optical states in the pulsed domain. In the continuous variable regime, 
the characterization of pulsed states of light is usually provided by time-resolved homodyne 
tomography (HT) [1-3], which allows obtaining, for instance, the Wigner function of the 
state [4]. Though HT has proved its effectiveness on classical and quantum states further 
optimizations to achieve a proper spatio-temporal mode matching of signal and local oscil- 
lator [3] are needed in the case of pulsed fields. On the basis of the results derived in [5], it 
was demonstrated that also direct detection can be used as an alternative to HT to obtain 
the Wigner function starting from photon-number distributions [6,7]. The only experimen- 
tal work [8] was performed with continuous- wave weak fields and single-photon avalanche 
detectors. The method consists in detecting the light exiting a beam-splitter that mixes the 
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signal with a coherent probe field, whose amplitude and phase can be varied (see Fig. 1). 
The Wigner function of the state can be written as 

r> OO 

W(a) = -J2(-l) n Pn, a , (1) 

*n=0 

where p n ^ a is the photon-number distribution of the field displaced by a. To obtain p na recon- 
struction/inversion methods must be implemented [9] that require a known and sufficiently 
high photon-detection efficiency [10]. However, in some situations, knowing the statistics of 
detected photons can be sufficient to characterize the signal field state, as demonstrated 
in the pioneering work of Arecchi [11]. In fact, if we are able to characterize the detector 
response, we can infer the photon-number distribution from that of detected photons with- 
out inversion. Obviously, obtaining or not the statistics of detected photons depends on the 
photon-counting capability of the detector, a fact that limits the range of intensities that 
can be investigated. 

In this Letter we exploit the linear response of a photodetector endowed with high gain in the 
scheme in Fig. 1 to directly reconstruct the Wigner function of some pulsed classical states of 
light in the mesoscopic domain. The detector output is analyzed through a new method [12] 
that allows obtaining the detected photon distribution without any a priori knowledge of 
the states [13]. We model the detection process as a Bernoullian convolution [14,15] and 
the overall amplification/conversion process through a very precise factor, which is taken as 
constant, 7. Our method has the advantage of being self-consistent as the value of 7 is ob- 
tained from measurements on the very field under investigation [12]. The idea is to measure 
a field state at different values of the overall detection efficiency of the apparatus, 77, and 
to calculate the Fano factor, F v = o~ 2 (v)/v, of the output voltages, v. Since we can write 
F v = (Q/n)v + , y (Q is the Mandel parameter), 7 can be obtained by fitting the experimental 
F v versus v data as all the dependence on the field is in the slope, Q/n. Once 7 is evaluated, 
we find the detected-photon distribution by dividing v by 7 and re-binning in unitary bins. 
Let us call p m ^ the distribution of detected photons, with (3 = y/rja the detected amplitude 
of the displacement field. If p m ,p is linked to p n ^ a by a Bernoullian convolution, we can write 
an expression analogous to Eq. (1), W{(3) = 2/tt ^* =0 (- l) m p m ^, where 

W{(3) = J d 2 p>e-^-^W(P'/Vv) ■ (2) 

is the Wigner function in the presence of losses [7] and W(P/y/rj) = W(a) is that of the 
photons. It can be demonstrated that, for classical states such as non-squeezed gaussian 
states and their linear superpositions, we have W(@) = W(a), i.e. the integral in Eq. (2) 
preserves the functional form of the Wigner distribution. 

In the scheme in Fig. 1, a frequency-doubled Q-switched Nd:YAG laser at 15 kHz repetition 
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rate (Quanta System) provides linearly polarized pulses at 532 nm of ~ 200 ns duration. 
The beam is spatially filtered and split into two parts serving as signal and probe fields. The 
photon-number statistics of the probe is left Poissonian while that of the signal is modified in 
order to get the different states to be measured. Signal and probe are then mixed at a cube 
beam-splitter and a portion of the mixed field is delivered to a hybrid photodiode module 
(HPD, H8236-40 with maximum quantum efficiency ~ 0.4 at 550 nm, Hamamatsu) through 
a divergent lens and a multimode optical fiber (600 /im core-diameter). HPD has limited 
photon-number resolving capability but is linear over a wide intensity range. Intensity and 
phase of the probe field are modified by a variable neutral density filter and by a piezoelectric 
movement, respectively. 

We describe in detail the reconstruction of the Wigner function of the vacuum state, which 
was obtained by simply blocking the signal. As the state is phase- independent, its Wigner 
function is rotationally invariant about the origin of the phase space and the reconstruction of 
a section is enough for full characterization: we thus only varied the probe amplitude \ f3\. We 
implemented the characterization procedure described above at each value of by changing 
r\ by means of a polarizer set in front of the optical fiber (we estimate a maximum value of 77 ~ 
0.31). We recorded N = 30000 laser shots for each //-value. The Fano factor for this series of 
data is shown in the Inset of Fig. 2(a): as expected for Poissonian light, the angular coefficient 
is virtually zero. Examples of the obtained p m ,i3 are shown in Fig. 2(a) as bars. The agreement 
with the theoretical Poissonian distributions (see Fig. 2(a), symbols+line), calculated with 
the measured mean values, can be checked by evaluating the fidelity f p = Y,m=oPmPmi i n 
which M is the number of elements of the distribution experimentally reconstructed. For 
the data shown in Fig. 2(a), we have / > 0.9999. By using the measured we obtain the 
experimental W{(3) displayed in Fig. 2(b) as dots. The experimental data perfectly lie on a 
section of the theoretical surface as calculated from the general expression for a single-mode 
non-squeezed Gaussian state [16], 

W{fi) = exp(-^^) , (3) 

7r2m th + l PV 2m th + l ; V ; 

in the case of no thermal photons (m t h = 0) and no coherent displacement (/3q = 0). For 
a quantitative estimation of the quality of the reconstruction, we evaluate the mean error 
e = Ek=iW th (f3 k ) -W sp ((3 k )}/N. The presented data yield e = -8.4 x 10" 4 . 
To enforce the usefulness of our method, we consider two further cases in which phase- 
independent signal fields enter the BS: (a) a phase-averaged coherent state, obtained by 
randomizing the relative phase between signal and probe with a piezoelectric movement 
operated at a frequency of ~ 300 Hz and covering 1.2 /im span, and (b) single- mode thermal 
state, obtained by inserting a rotating ground glass plate into the pathway of the signal 
field, followed by a pin-hole selecting a single speckle. In the main panels of Fig. 3 we 
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present the p m $ (bars) obtained from two sets of data corresponding to different \/3\ values. 
The reconstructed p m ,p for selected values of \(3\ were then used to recover a section of the 
Wigner function (dots in the Insets of Fig. 3). 

To compare the experimental results with theory, we have to take into account the mode- 
mismatching between probe and signal. In fact, the imperfect superposition at BS leaves 
some residual probe light that is measured together with the mixed state. Thus the measured 
distributions result to be the convolution of those of the mixed state with the Poissonian 
distribution of the residual probe, while the Wigner function is the product of the expected 
Wigner function with that (Gaussian) of the residual probe [8]. In ) the overlap 

between signal and probe was ~ 91% while in case (b) it was ~ 65%, owing to the further 
difficulty of producing a stable single-mode thermal state. 

In Fig. 3 we plot two theoretical curves superimposed to each set of experimental p m> p. One 
(white symbols+line) was evaluated from the expected theoretical distribution in which the 
measured mean values for signal and probe were inserted. The other (black symbols+line) 
was evaluated by calculating the expected theoretical distribution with a probe reduced by 
a factor corresponding to the overlap and then convolving it with a Poissonian distribution 
having a mean value equal to that of residual probe. The latter distributions exhibit better 
superimposition, also testified by the increase in the fidelity from / ~ 0.998 to / ~ 0.999 
in (a) and from / ~ 0.990 to ~ 0.999 in (b). Finally, in the Insets of Fig. 3 we plot the 
theoretical Wigner functions obtained upon the correction procedure described above. The 
non-corrected expected surfaces are not plotted for the sake of clarity. The expected Wigner 
function of the phase-averaged coherent state is 

W{(3) = -JoWH&l) exp[-2(|/3| 2 + |/5 | 2 )] , (4) 

7T 

in which Iq(x) is the modified Bessel function and \/3o\ 2 is the mean number of detected 
photons in the signal at the output of the BS. The Wigner function of the single-mode ther- 
mal state comes from Eq. (3) for /3 — and m t h equal to the number of detected thermal 
photons in the signal field at the output of the BS. Since during the experiments we had 
some fluctuations of the signal field intensity, we kept the value of the signal as a fitting 
parameter and compared it with the average of the measured values. For the two cases un- 
der investigation we got \Po\ 2 = 1.41 and m t h = 1-96 to be compared with 1.55 and 1.63, 
respectively. The errors on the Wigner were e = —7.8 x 10 -5 and e = 4.1 x 10 . 
In conclusion, we implemented the reconstruction of the Wigner function of some pulsed 
classical light fields starting from direct measurements of the statistics of the detected pho- 
tons. What made the experiments feasible was the choice of a single-photon sensitive linear 
detector together with a new data analysis method. Works are in progress at our Labora- 
tory in Como, whose aim is the reconstruction of a phase-dependent signal field, namely a 
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coherent field, as well as to characterize non-classical states, both squeezed and conditional, 
in the same intensity regime we explored for classical ones. 

The Authors thank M. G. A. Paris (Milano University) for fruitful discussions and P. 
Salvadeo and A. Agliati (Quanta System, Italy) for the long-term loan of the laser and for 
the promptness of their technical assistance. 
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List of Figure Captions 

Fig. 1. (Color online) Scheme of the experimental setup. HPD, hybrid photo-detector; BS, 
beam splitter; NF, neutral density filter; P, polarizer. Components in dotted boxes are alter- 
nately activated to produce different signal states: V, vacuum state; T, single-mode thermal 
state; AV, phase-averaged coherent state. 




Fig. 1. (Color online) Scheme of the experimental setup. HPD, hybrid photo- 
detector; BS, beam splitter; NF, neutral density filter; P, polarizer. Compo- 
nents in dotted boxes are alternately activated to produce different signal 
states: V, vacuum state; T, single-mode thermal state; AV, phase-averaged 
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Fig. 2. (Color online) Vacuum state, (a) Inset: Fano factor, F v , as a function of v. Main: 
reconstructed p m .j3 distributions at different |/3|-values (bars) and theoretical curves (sym- 
bols+lines); (b) experimental reconstruction of a section of the Wigner function (dots) su- 
perimposed to the theoretical surface. 




Fig. 2. (Color online) Vacuum state, (a) Inset: Fano factor, F v , as a function of 
v. Main: reconstructed p m ^ distributions at different values (bars) and the- 
oretical curves (symbols+lines); (b) experimental reconstruction of a section 
of the Wigner function (dots) superimposed to the theoretical surface. 
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Fig. 3. (Color online) (a) Phase-averaged coherent state. Main: reconstructed p m ,p distribu- 
tions for two values of \(3\ (bars) and theoretical curves corrected for the overlap parameter 
(black symbols + line) or not (white symbols + line). Inset: experimental reconstruction of 
a section of the Wigner function (dots) superimposed to the theoretical surface; (b) same as 
(a) for the thermal state. 




Fig. 3. (Color online) (a) Phase-averaged coherent state. Main: reconstructed 
p m; p distributions for two values of (bars) and theoretical curves corrected 
for the overlap parameter (black symbols + line) or not (white symbols + 
line). Inset: experimental reconstruction of a section of the Wigner function 
(dots) superimposed to the theoretical surface; (b) same as (a) for the thermal 
state. 
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